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^ It  Is  (ie.‘3ii'S'blc  to  predict  scsa-curface  tcEpeistuies  into  the- 
Predict  loa  iatenrale  HSBy  1)3  a fev  dsys^  a fsi'  acK.nthfc;,  or  even  ceveral  years 
for  -loag  terse  prediction » Hiero  exist  several  :?.ong  tera  recerde  of  daily 
surface  teo^peraturea  taken  at  fixed  IccationSc  lijese  records  veay  fx-esm 
five  to  ei^t  years  in  length  with  only  a awderata  araouat  of  njlseijog  data, 

Ube  ultiuate  goal  of  a current  study  ie  the  ceteiiair5i.'f  loa  of  a method  of 
using  certain  characteristics  of  these  long  tm-Ei  reccids  txj  forfccawsl,  surface 
te^^pej-atures  several  moaths  into  the  future^  The  seine  aathods  could  he  ex- 
tended to  forecast  tactieeJly  iE^portaat  undeniatsr  sound  velccitles  as  a 
function  of  depth  and  geognphical  location <,  iortsdiate  goal  of  this 
sjesaorandufli  is  to  record  e hackgrouad  for  such  & study,  to  do  cert«.;ln  pre- 
ilmlnarj'  work  in  tliae  series  aaelysis,  end  to  outline  tlje  neact  phase  of  the 
worko  N 

A BACJaiiiOTKD 

Pre'/ious  ivork  at  NEL  by  E«H.>  Aiiderson  and  3.  Wo  Snedecor  has  included 
inultlpie  regression  analyses  of  surface  sound  velocities  in  the  region 
30°  - B,  120°  - 123*^  W for  the  fiscal  years  1951  to  19550  In  parti- 
cular, regression  equatltMis  as  cca^lex  es 

where  D ie  time,  G is  loagitude,  T is  latitude  and  V is  surface  eoustd 
velocity,  have  been  fitted  to  observations  for  the  Individ'jal  years » Tbe 
above  regression  e<iuation  considers  the  mtin  effects  of  tiiae,  Icngitude  and 
latitude,  and  the  interactions  of  tliae  with  latitude  and  longitude  with 
letitudeo  It  should  be  noted  ths,t  a fifth  degree  polynoinial  fens  been  used 
to  rerpresent  the  variations  vdth  respect  to  time  during  a year^  Eiis  approach 
differs  from  that  of  the  present  jaeasonandun  in  which  periodic  functions  of 
tlaie  are  usedo 

Oceanographic  data  for  Oceeui  Weather  Station  PAPA  {50°H,  li^5®W)  have 
been  analy&ed  r&ther  extensively  as  pert  of  the  Environiisaital  Statistics 
Stodics  (I!FE;  Problera  Zl-Sl)  jho  d«to  aro  beinn  oolloctod  by  Conaclians 
of  the  Pacific  Ocesnogrciphic  Grr.up,  ISauaiao,  B,C  , and  axo  awailablo  in  A 
seq.usnce  of  reports  published  by  the  Pioberics  Peaearch  Board  of  Canada 
In  their  Manuscript  Report  Series  (Oceanographic  end  Ltiaiologlc&l) « The 


dilva  incdUiie  periodic  s’eadinos  o^r  Jvith'i'theMaograph  ftixsts,  surface  ten^xsra- 
ture  laeasuxeaentsj  and  surfece  and  depth  roeaEuroments  of  salinity,  dissolved 
oxygen  and  silicates^  IWice  daily  bathythexioograph  casts  are  made  at  0200 
and  1700  and  daily  sa&iples  for  surface  saxinity  determination  are 

collected  at  0200  G»MoT.,  This  maraorandtJSQ  Is  concerned  with  the  long  term 
time  series  analysis  of  the  daily  surface  ten?»erature  data  as  obtained  from 
hathuhersnograph  record  So  The  0200  data  are  U8cx]  vihpo  avallablo}  if 

not  available,  other  data  for  a given  day  axe  usedo 

In  addition  to  the  Pacific  Ocean  Station  PAPA  data,  which  received  the 
initial  intesnsive  analysis,  there  exist  several  years  of  data  for  Atlantic 
Ocean  Stations  DELTA  (440Ns  4low)  and  EaX)  (35®N,  480W)  which  have  been 
analyzed  in  the  same  manner  as  that  for  PAPA*  Tine  PAPA  data  have  been  collected 
by  oceanographers  for  specific  oceanographic  uses,  and  are  definitely  of  higher 
quality  tlian  the  DELTA  or  ECHO  data^  The  Station  PAPA  data  are  for  the  time 
interval  1 July  1956  - 30  Juno  1961,  Station  DELTA  for  1 September  1949  - 
31  August  1957,  and  Station  EOK)  fox  1 September  1949  - 31  August  1956o  The 
DELTA  and  ECHO  data  were  obtained  from  the  Hydrographic  Officeo 

TIME  SERIES  DEOMiPOS  ITION 

Any  record  of  surface  temperature  necessarily  reflects  soma  sort  of  an 
annual  oscillation  and  possibly  some  sort  of  a trend.  It  is  convenient  to 
rmnove  the  annual  oscillation  and  tlie  tr^xl  from  a series  of  observations 
before  applying  any  statistical  prediction  technique.  The  resulting  surface 
temperature  anomaly  becomes  the  dependent  variable  for  a time  series  analysis. 

If  the  anomaly  oidiiblts  a useable  persistence,  then  a con^-osite  prediction 
formula  consists  of  three  terms, 

1,  an  estimate  based  on  the  trend, 

2,  an  estimate  based  on  annual  oscillation, 

3,  a statistical  predictor  based  on  the  anomaly  characteristics. 

The  general  concept  of  trend  is  that  of  a gradual  change  in  a system  over 
a long  period  of  time,  where  "long"  is  a relative  term.  An  essential  part 
of  the  concept  is  that  the  change  In  the  systms  over  fairly  long  periods 
is  monotonlc.  It  is  always  possible  that  a trend  can  be  confused  with  a 
slow  oscillation}  this  is  certainly  true  if  the  record  length  is  short 
compared  to  the  period  of  oscillation,  If/lth  respect  to  surface  tenpera- 
ture  prediction  a few  months  into  the  future,  such  confusion  is  unixqportant. 
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For  Ihii  tirj2  muier  sor.‘vicl)?.rr.ti<xe  .tii  il-i'i.a  :jsa 'o;fa,'id-.ifli,  tLc  rccci-a 

lacfttbs  of  fiv€i  -CO  eS^'it.  ystass?  SMSf  ra/ohar  1’’  Hii>s  -^or  arjr  tr-aiKl 

deterrsiaaecidi,  parb?.e?,i;!iu;:i.y  in  Hght  of  tav  cxist-iintcij,  to  hi  aotn.  ir-.tor.,  of 
a very  doiaiiinct,  SAjauft.!  os*:i:Llat'-<sx-.  ivi^iLou-i;.  galas  cixj:'fca.il,  cay  diaasgec. 
in  the  systsia  s«cm  to  be  yea?.'  to  year  r-^^D-utt-as-tlana  not  the  Jtoag 

tens  req.uirefjc»itB  of  treiido  Fox  the  foi:loKin3  dlss-nxfssioi;.^  it  1®  risEixcsed 
th&t  xjo  tx-ead  exi&tB  in.  the  e’xri’t.cs  t'Ksjpera.tux'O.  raeordso 

m>L£sm  CF  THi:  mr/os  safa  j m.mAh  of  m:  Ammi  vArctAnoii 

Tee  Stetioa  FAFA  d.ote  eoaSyzed  ia  this  5''.MaorS!i’3duc\  cenaxtst  of  .'.,136  of  tha 
poeolble  I826  Burface  teii33erc.ta’’e  xaeaaiu'essanta  for  the  ttao  Interrval  j Tuly 
1956  “ 30  stoe  195io  Except  for  tixo  j^ps  of  about  cxiS  r^atli  the  .aisalrtA 

data  dayc  ers  di'.jtribute<i  laore  ox''  less  svt  roudoai  wltbia  the  five  ye&i'  peviQd.„ 
?S.gure  1 display®  the  FAPA  date,  •aodlfled  f03.'  plotting  convEHience  to  the  ex- 
tent that  the  isedlBn  Ci5aez".>al-,icio  in  ej-ch  successive  five  «iay  lUiterve.!  has 
been  plcr}’.tedj,  rather  thffjs  the  ind:f.vldyal  c^AHervP-tioaSc  It  sfceroiii  be  eephS' 
sised  tljat  ell  aBaJ.ysic  was  actual3y  psrforaed  cn  -khs  individiial  daily  cb- 
eerviitlaiSj,  aad  not  oa  these  plotted  Jtsed3.ai;i8<. 

•Hie  first  dtevicni®  feature  of  the  PAFA.  chxoa  S.e  the  ennur.!  oscD.Iatton 
xfitb  fairly  emirp  aaxiiKi  in  A'jgust  to  Oatdfcc**'  eud  rather  flat  0i3i.h2R  in 
Fobx'uax'j-  to  lAx;',  Sines  tlvs  dcata  arc  ©esillfftory  in  nature  with  a poriod 
nesessfirily  close  to  one  cslcudar  year,  it  ia  JuBtifi&hle,  iVesra  a ctatistical 
point  of  viert?,  to  fit  to  t5ie  diate  sck;?  theoi'ctical  fu;ction  which  oociilatss 
with  period  cne  yesxo  Such  a fuiioticx-  wcuid  be  fitted  xo  tlic-  d&ta  either  to 
pro-.'ide  an  xmderstsadiag  of  'ihe  oscil-latoxT  nature  of  surface  tci%)&rc.iuraB, 
or  to  provide  a prediction  foiraiLs.  for  there  toiigKiiaturec^  Ihe  siiaplest 
such  fimctiou  l«j 


/' 


•--■  --  ;Si  5 > 'S:..C-.OS'J&--  O , 


jt-yo 


(lA) 

(IB) 

where  D is  tlas  aeanured  in  days  fixo  e.rbitrcxy  or.iglnj  <jz  the  fitted, 
value  of  ourl^e  tcaper^turej,  and  either  ^ i>'.  or  e(iuiya3.eT;.v.lyj; 

l^j.j  e,:i.*e  refcrcssiou  coe.rficienta  to  be  estiriete-d.,  ®ie  ^rear's  length,  hore 
shovn  ao  3o5  dii;;^^.  '■•Jouid  be  wrltt.ju  to  what,f-/er  sy'cniracy  is  ,1uct?.fi«d  in 
the  ertur.i  co;;53irt*t.ians  0 

A convantiono-l  way  of  fitting  such  a £\a.ctlon  B'p-ntica  (3.)  to  the 


3 


joints  of  if'iS'-'-va  1 is  tiic  &‘-rchO'i  of  L-jjist  SQVU?.rai5 « It  la  dectcledly  esMrler 
f,o  ues  IS<^uatio:i  (IB)  th»n  (l/O  ii?.  tiiia  aiafchod.  of  fitting  sisiee  ^cj 

j ~:i.  wiLl  !Xiipfic.T  as  tbe  solution  of  three  siasi'jtaneouc  .linecu'  ©(luatiooe, 
•hjereas  ■t^jould  atipeai'  be  tbe  solution  of  three  quite  ccss^pllcftted 

xijsi^ltaneous  £-csi"liKear  viquE.tiouso  Ivbz-a  specifically,  the  Bethcd  of  Least 
Sqixsree  consists  of  BsJ.uljais'.iag  with  respect  to  the  subs  of 

squares  of  'the  dlffareacec  of  ctssrved  atid  predicted  values 

Q - 

Z (K—  3c  i 

'fhsre  the  BUiaavtion  is  over  tize  E seEiple  poiJits^  Thue  the  equis.iirsja 

■§' Jr.  -•*  O i -k>r  k 0, 1 I 

are  to  "be  solved  BxssuJ.tancou83.y<,  After  sass  Msanipulatioa  these  eoiaitioeui 
ca®  he  Witten  a® 


Ptj 


+/Si2,s;rig  D 


D “ 2-(X 


(2) 


Btviog  determined  B, , 

-•  f 


it  is  easy  to  ssLtata 


!.?r 


CJJrc^Td'h 


jSE. 


(3) 


Before  d;lflcusslng  the  numerical  i-sBultB  of  fitting  the  PAPA  datft,  cco- 
Gldei  the  folloif:ir*gc  Althougjh  s viEnj®!  ofoservatton  of  figure  1 vould 
suggest  the  fitting  of  Eenjatien  (iB)  to  the  data,  and  is  a rreaBonable 
technique  in  an  Investigative  problem,  it  is  perhaps  enli^tenljog  to  present 
the  PAPA  date,  la  terms  of  tbe  autocorrelatioa  fuactlcn 


R4  - 


cx'V  X\ 


— ... 

^fiia 


O/’  1 • ^ J 


where  all  the  suKeo-tiotn  sigas  Mean 


(iv— i 
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k~  I 
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LAG  j (DAYS) 

AUTOCORRELATION  FUNCTION  Rj  FOR 
THE  ORIGINAL  STATION  PAPA  DATA. 


FIGURE  2 


L 


i 


She  aut^ori\?l^l:io3a  "urxewrxa  o^^fiiiSd  tia^  R i lor  Jiix 

o ^ 

!iierO;,  nafiuljg  tbRt  the  tjUoe  serica  ii3  perfect iy  coxrelxcQd  witL  itt^elf. 

If  the  PAPA  tlae  series  repeated  itseLf  sxactly  OTery  ./eary  t>3a  ca^tocorre'tp.- 
tion  for  :j  « 1*  2,  years  tfcsiM  a.i^  be  = 1,  siacs  there  vould  be 
perfect  correlation  wliesnever  tie  tine  3crie3  arjd  tie  iasged  tfjos  series 
were  exactly  in  pbaseo  Since  the  dat'-i  do  not  x-epeat  theKseives  cx&ctly 
every  yesr^  the  autocorrelatioa  fyucctioa  csiJ  never  be  erne  except  at  j » 0, 
but  as  is  seen  frea  Figure  2,  tie  autocorrelXwticKi  functloei  las  very  large 


pooitive  peaks  at  yearly  iaterva:LBt  ‘ais  is  a defrLults  iiidic;4t}.on  tbP.t 


there  are  oscillaticn®  with  period  osxe  year  in  the  PAPA  detsi.. 

Pitting  Kquatioc  (UB)  to  the  data  yi-eldjE  the  regression  equation 

A "'  3. ^'•3  O — 

G.54-  '‘'.3.8-“^ 


A 

where  the  units  of  X are  degrees  Centigrfiuioj,  suid  day  I>=1  ia  Janixaiy  1» 
As  measrures  of  the  goodness  of  fit  of  this  eciiatio«i  to  th®  actual 
data  we  can  consider  three  related  qufvntltiea * 

Ic  ffiia  Bazltiple  corrslaticsa  coefficieaat  between  the  obssrveo.  sad 
fitted  X’So 

2«  llhe  variance  of  the  obsei'^/a.tijias  e.bw.t  the  fitted  curve. 


(5) 


3.  The  P»ratio  of  the  regression  sias  of  squaree  to  the  uaej^lained 


BU£ti  of  squares  o 

Cocsidering  these  quantities  in  turn,  the  equ&r®  of  the  SKilt?.ple 
ccrrelatlca  coefficient  bett^ea  the  observed  surfece  tis^jsrature  X end 
the  fitted  surface  teB5>eraturc  for  the  PAPA  data  is  = 0.863#  freas  which 
Rt>0.929.  The  quantity  R is  a sseewsure  of  the  strength  of  the  linear  relation- 
ship between  the  two  variables.  A value  of  R as  large  as  0.^9  J3s®>li08  that 
Equatlaa  (5)  provides  quite  a good  fit  to  the  PAPA  data. 

ia;o  variance  of  the  observatioce  aboirt  thoir  mean  val'ie  is  8 504.  The 

variance  of  the  cbser/atlons  ffihOijit  the  fitted  cui*ve  is  1.220.  The  reiatfcaa- 

o 

ship  between  these  ■'/arisneen  Ciid  ir  ic  obr;c-ive<l  to  be 


jj^  ^ 30 


•fehioh  leads  one  to  osodsoea 


0. 8 5-3  j 

.xs?  re  being  the  of  vari^ce 
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STATION  PAPA  DEVIATIONS  AFTER  ANNUAL 
OSCILLATION  HAS  BEEN  REMOVED. 


FIGURE  3 


FISCAL  YEAR 


exp^lned  by  fitting  the  j'cy;pessi(Ki  cu3r/e«  Tha  Btendard  deviation  of  the 
obeBrv’ations  about  the  regi'caalon  cuwe  1*  ‘'/[.2^?,0  *'•  L!OS**C.- 

Pinaliy,  the  total  e\m  of  Bqm.ceB  of  the  PAPA  cbsei'vationa  about  their 
B-'san  xB  10^114,6.  ‘£he  reduction  in  sum  of  aciuares  due  to  fitting  Equation 
(5)  is  8729*2,  vhlch  leaves  en  unexplained  Bvja  of  squai-ea  of  I385.4,  With 
the  proper  eorpjftllty  eflsusptioao,  the  last  tvt>  euaia  of  squares  are  Independently 
Chl“3quare  distributed  ’./ith  2 and  II33  degrees  of  freedom,  respective?  y,  so 


that  the  F-ratio 
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~ rrs- 
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cin  be  vrlttea  end  interpreted*  !Shis  value  of  F is  highly  significsuit 
since  the  1^6  critical  value  of  F for  the  above  degrees  of  freedora  is  FC.Ol  ■ 
4*6a»  One  concludes  that  the  regression  Eque-tion  (5)  significantly  reduces 
the  total  suffi  of  squares*  The  analysis  to  foij,ow  will  shas  that  the  \uiex- 
pltlned  sua  of  squares  can  be  reduced  f\utker  In  a sseaningful  way. 

The  clue  that  an  annual  oseilletloa  existed  la  the  PAPA  data  was  pro- 
vided by  visual  6b3eiv».tlon  of  the  original  data,  or  of  the  autocoi-relation 
function.  Statisti<:»l  Justification  for  fitting  a given  functiem  to  data  is 
provided  a-posteriori  by  the  F-mtlo  saalyels,  or  its  equivalent. 

AEALYS.TS  OF  TIE  SSATIOE  PAPA  IlATAj 

WimAL  OF  THE  SHJflASttfJAL  TARTATIOH. 

'Eae  fit  of  Equation  (5)  to  the  PAPA  data  can  perhaps  best  be  displayed 
visually  by  plotting  the  deviations  of  the  obasr^.-atiens  from  the  fitted  curve, 
an  he»  been  done  in  Figure  3*  This  persalts  an  e^»anded  vertices,  scale  over 
that  of  Figure  1.  Again  the  teedi^  observation  in  each  successive  five  day- 
interval  has  beesn  plotted.  A feature  of  the  deviation  plots  is  the  existence 
for  each  fiscal  year  of  tw'o  fairly  well  defined  msxiaa  at  rou^ilj.'  cne-half 
year  intervals  separated  by  fairly  well  defined  miniaa.  For  i960  the  first 
jnsxijHum  occui«  in  the  last  5 days  of  1959>  Figure  3 provides  visual  evidence 
of  the  existence  of  sane  semiannual  osciUatior  in  -the  PAPA  data, 
as  fin  oltoznoto  due  to  tho  mdatonco  of  tho  ssdannual  vasnation;  00 
again  zsnke  use  of  the  autocorrelatioa  function  as  was  done  in  Figure  2. 

Bie  autocorrelAtlon  function  of  the  devla  tlons  of  Figure  3 is  casqn-ted 
end  displayed  in  Figure  k.  This  figure  showr  quite  clearly  that  a semi- 
annual. oscillation  exists  in  the  PAPA  da-ta,  althou^  -Ute  peaks  in  Figare 
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FIGURE  4 


h are  «i?aiJLor  in  Esaig'3.ltuc',e  cisci  tii?.or,c  of  Figazra 

HaviBg  diSBOnstmtcd  iise  !aiis5i;«xce  ox  ooEje  sort  of  a tcaidsnual  cacilla- 
tloa  ixi  tiie  ds.ta,  it  is  nw  co.,ive»i.aiit  to  revarb  tc  th&  original  PAPA  S/ita 
and  reaove  the  r-aznial  »nd  soa5jSBr.y.el  oseiliiitloru;  fflEailten^cualy.  Analogous 
to  tlK  woii  above,  the  3:iaEileet  function  cou'^irilng;  both  enaual  end  semi- 
annual  o8Sl3Jls.tor:y  teroti  ie 


Ssn^^lD-e,)  -b  D~Qxx)) 

= ■'.'7^1  s:iT^^  o A-?. 

^ c*  -4-  ^ D . 

ffiie  qxjcntlty  Q o “ii  (X-X^)*"  Is  vlth  I'espect  to  the  five  « 

•Jtot  ie,  the  aifrations 

ac?  , - 

apX  " 3 3 2.J  3^  4^ 

are  solved  sSaalatiieoueiy.  Analogous  to  SjuationE  (2),  there  ere  five  linear 
equations  in  the  to  be  solved  siEiultsaeoizsiy. 

The  Hdalmuai  of  Q obtained  by  siiail.tfmeously  fitting  tcith  anauai  and 
3«aiian3ual  terjaa  is  the  ijajce  3aia3j5!um  that  iroald  be  obtained  if  a fit  of 
aenu&l  teres  were  followed  by  a cepaj»te  fit  of  sa^rLannual  t-arrs  tc  the 
renultiug  deviatiens,  obvious  c2itaig3S  of  origin  be:jt>g  EBcleo  ‘Bie  two 
separate  regrecaion  equations  added  toge'cher  vrould  be  tbo  saar  as  the 
sisulteneous  resressi.0Q  equ&ticco 

B»e  eaplitudee  and  phacas  are  given  by 


m 

(6B) 


C<,  2^ 


(7) 


Fitting  Equation  (6e)  to  the  dats^  yields 


/. 


'P  4'-  - 


3,3  Sin-I&'D 


i.c^2  0 


'‘a£>s_  3'os^ 

-f-  on  b-O  n P -■  UA05  cos^  D , 

« a4<5!">’aa^,  ^i-^~^^(0-[oZ.0)--r0.33 


(8) 
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vhere  X i«s  5»  CectigWide  esd  )>«1  Is  Ssic:uiay  1* 

The  3q.\iiare  of  the  sultlplo  correlAticn  eoefflcleot  hetweffis  the  dbeerved 
eurltice  tasg)er«.tiire  X ecd  this  nev  flttcfl’.  X^  is  S^aO.903f  ud  vr 

jjccreeoe  fraa  £^«0.863,  aod  1^0. 969,  vhsn  caUy  the  annite,!  oeeillAtioc  wu 
fitted 4 

The  vArlwoec  of  the  ohservutions  saMut  the  fitted  eurre  is  0.865»  e 
decrease  .•?2‘ca  1.220.  The  eerrespoBding  staxidard  denristlan  Is  0.930^C,  a 
decrease  i'rcji  1. 105^0. 

More  the  addltloael  reductlcn  Is  sxxe  of  squares  due  to 

fittisfi  the  8(!al«cBV!al  oscilXatioa  is  403.0,  lesvlng  ea  uncaqplained  ana 
of  squares  of  9^.4,  aad  resulting  in  on  F>ratlo  of  232,  still  hl^ily 
ai^ilflceat  at  the  1$  level.  At  this  poi>at  it  vould  sees  aj^ropnate  to 
display  a table  alxsllar  to  a&  AnalyBis  of  Vhrlaxiee  Table  in  order  to  teat 
the  sigoificaoce  of  the  reductica  in  sthc  of  squares  due  to  each  additional 
regression  varlshle.  liRiis  is  done  in  Table  1< 


TABLE  1.  VPiPk  ASALXEIS 


Variable 

Sum  of 
Squares 

]Degreee  of 
Freedom 

Mean 

Square 

P 

Ratio 

Total 

10114.6 

1335 

Annual 

8729.2 

2 

4365. 

3569 

Une:q>lAised 

1385.4 

1133 

1.223 

Sesiianaual 

403.0 

2 

201.5 

232 

unexplained 

982.4 

1131 

0.069 

(Significance  Levels;  F0.05  » 3*00,  FO.Ol  * 4,62) 


On  the  basis  of  the  above  table  one  concludes  that  the  regression 
Equation  (8)  reduces  the  total  sun  of  squares  slgalflcantly  aore  than 

that  acsoual  aad  aeoiannual  oscillations  should,  be 
data  alaultaneously  prior  to  further  analysis. 
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does  Equation  (5),  so 
ranved  trosi  the  VAPA 
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STATION  PAPA  DEVIATIONS  AFTER  ANNUAL  AND 
SEMIANNUAL  OSCILLATIONS  HAVE  BEEN  REMOVED 


FIGURE  5 
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As  before,  the  fit  of  BquBtioo  (8)  to  tho  PfiPA  date  cmi  be  diBpJjtyed 
visually  by  plo-ttiog  the  deviAtiOBs  of  the  observaticr.E  frosc.  the  fitted 
curve,  as  has  been  done  in  Figure  5»  Alternatively,  the  scutocorrelatloa 
function  of  these  deviatioas  has  been  plotted  in  Flg'Are  6»  Heither  Figui*©^: 

5 nor  6 indicates  the  presence  of  any  osctUaticn  Tfith  pericd  loss  tbna  ont 
year  in  the  deviations  after  the  annual  and  eerslvnaual  CBclllsticna  have  been 
reasoved.  ,!f6  further  analysis  of  the  unexplained  s^^aa  of  squares  wiL'.  bo  attaai^ted 
In  this  nwaoiandias,  Ibture  aoalyBis  of  the  rssiduel  variability  vill  be  outlined 
in  the  fiiatl  secticBi  of  this  nteasorandura^ 

AmisBxs  OF  8!2»£nt>a  mm.  asd  sisatsou  bcqo  data 

file  aethoda  developed  above  for  the  aatalysiR  of  the  FA?A  data  have  beai 
applied  to  the  SSlfiA  a^d  ECK)  data.  Similar  resulte  bHve  be&ti  obtalced.  fiie 
reoultc  for  the  three  statlouB  ca-e  preseated  below  in  trhulsr  form.  Table  2 
preseatG  the  regression  coefficients,  the  iaultiple  correlation  c-oefficiants, 
and  the  standard  deviaticos  of  the  oboervatic«is  about  the  regression  curves. 
Thbles  3 find  4 axe  Analysis  of  FarlSKice  tables  correepondiag  to  Sable  1 for 
the  PAPA  data. 


sasLs  3.  mm  Amssm 


variable 

Sua  of 
SQ,uares 

Degrees  of 
Fraedon 

1 .11 

j'teoia 

Square 

P 

l:htit.e.O 

Total 

9R85 

1287 

Anssial 

6203 

2 

3101 

1293 

IbexplaSned 

306E 

1285 

2.398 

Seadannual 

2 

285 

ll*6 

IMexplAined 

.1283 

1.958 

TABLE 


Units  of  regression  coefficients,  amplitudes  and  standard  deviations  are  degrees 
Centigrade.  Units  of  phases  are  days. 


4.  Bern  misjsxs 


l^irieble 

Sun  of 

Degrees  of 

Mean 

F 

Squares 

Arecdoa 

Square 

Batlo 

iJotal 

13625 

1530 

Annual 

12065 

2 

603^ 

5902 

Uherplained 

1561 

1528 

1.C22 

Seaianmial 

493 

2 

2lt6 

351 

Oaezplalned 

1068 

1526 

.TOO 

(SigBlficance  leverus:  10-05  =>  3-00,  FO.Ol  • 4o62) 


OnSCLUSIOSS 

Ttie  ias>ort&st  dlscovex^  described  in  this  anostrenduai  is  the  hiifhly 
sl0iifle«]t  reduction  ia  mxaa  of  aqissres  in  surface  teaqperetures  attributable 
to  the  fittiag  to  the  data  of  anriunl  and  eesdannual  osciUatoory  terse  o Coe- 
aiderlng  5M>le3  X,  3 wd  4,  the  large  F>mtios  corresponding  to  fitting 
acadannual  tezass  even  after  the  highly  slgnlfieent  annual  tene  have  been 
fitted^  iadleetes  that  both  azasual  and  seniaonual  oaeiHatloBS  should  be  re- 
iaoved  froB  surface  temperature  data  before  statistical  analysis  of  tbs  aasnallea 
is  perforated. 

TUraing  actentlon  to  fable  2,  ahleh  ecoeems  the  eoaffielsnts  and  fit  of 
the  n^ression  curves,  a feature  of  the  table  is  the  slvllarity  for  the 
three  stations  of  terns  aithln  each  of  the  coIumib  azeept  ^ and 
fbese  stnilAritles  exist  even  though  tlM  stations  are  in  differmt  oceans  and 
even  thou^  there  are  dlfferenc^M  in  the  quality  of  the  recorded  data.  It 
would  appear  that  the  najor  differences  anong  stations  is  accounted  for  by 
the  differences  in  absolute  magnitudes  reflected  by  the  . fhe  dlffar” 
encas  nnoag  ststlans  of  the  angpUtudeB  are  relntively  ninor,  as  are  the 
differences  of  the  pfaasea,  althout^  it  la  the  differences  anong  the  pfaaaea 
Og  idilOb  detenBiiie  the  fine  details  of  the  shapes  of  the  individual  beat 
fitting  curves. 

Qoe  concludes  that  the  annual  variation  in  sea-surfsce  tsmperaturet  is 
wall  represented  by  tlie  sub  of  two  very  siaple  functions,  a sine  function 
with  period  one  year  sod  a sine  function  with  period  one-half  year.  Al- 
though an  annual  oselUatlaa  In  surface  temperatures  eas  certainly  aipeeted. 


10 


^ 1 

the  erUrteace  of  the  aea<lMacxiiil  m<:il3^±ioa  naa  cor,  erpec-i-ed,  | 

1 

\ OP  PtWOllB  WCSZ  I 

^Btls  aKstojao'lum  hao  '’oeoonstratcd  tlis  exi.Btesce  of  ajinv-C-J  and  aesilacaual 
oscilletions  in  aurfttca  tatroerature  data  for  oi’.e  Raclfic  Ocena  and  two  Atlantic  | 

. Ocean  statlcsis*  When  these  osciUatioae  are  recoved  froui  the  dAata  there  re-  ‘ 

»in  tlae  cerles  cossistlns  of  the  deviations  from  the  resreosioo  curves.  i 

j Bef erring  to  Figure  6,  which  nhoifs  the  eu.tocorrK'ijatlO!a  functioo  for  these  | 

derrlAtlons  for  the  PAPA  data^  there  are  indicf’.tiono  of  a perclsteice  of-  BOi'/,  j 

oef-nsoty-five  days  in  this  data.  ‘JJhe  w.luea  of  tmtoeorrolatim  cait  to  this  i 

I point  are  prcSsahly  elgnif leant,  'fhis  perslstfince  also  earaoars  in  Figure  5,  J 

j shoivlag  the  deviations  thenmelveo,  Iji  that  Buccessivs  deviatiens  do  no’s  vary  | 

Irsndxnly  about  the  but  ore  above  or  below  the  oxio  for  fairly  lOHg  i 

periods  of  tiiua.  It  is  fiuestionable  that  idiG  oscillatiocs  in  the  auto-  ! 

correlation  function  of  Figure  6 are  5jaj33iingful  jHuch  beyond  100  days.  They 
are  due  to  osore  than  the  characteristics  of  the  p9r'ticular  five  ! 

I year  sample  enrSULahle  for  analycis.  | 

^ If  a Persia bance  of  two  or  three  Bontha  exists  In  a time  series ^ then  ! 

I it  should  ba  possible  to  use  this  parsf^tsnee  to  predict  into  the  future.  I 

Present  plans  envisage  n prediction  fonaula  of  the  for.a 


Y 


' ^ n -!y-  i-J 


0(j,  <-'u  '''ij— !; 

where  Y is  the  surface  teD5>eratur3  snoaely.  tJhiD  equation  indicates  that  a 
teaflperature  anaaaJy  k days  in  the  future  will  be  predicted  as  a linear 
CQ(ift>ljQatlon  of  the  zn-l  most  recant  observsticais.  With  k end  n as  psrs^neters; 


the  Cj  will  be  determined  by  Kdnimlalns 


A />'!  , for  a B-ztown 


'ivlS  ''>.>!<- 

portion  of  a time  series.  Tba  ultimate  evaluation  of  the  predict!^  fcimula 
will  be  based  cm  hov  veil  It  pi^dicta  the  TsnknoMu  (to  the  predictor)  future  of 
the  tlaa  serloo.  As  mentianed  esrly  in  the  Baaac7eRd\m:i,  the  eoniposite  pre- 
diction foxiaila  would  caabizja  the  above  snco>aly  prediction  with  trend  and 


ItUittricsd^  ei^cai  u physical  ia-i.«rprtjtfi.v.ioifli  iS  difficult*  The  atatioca 
diflcusefed  iTi  this  ine!icoia«.dum  are  ia  tae  opcm  oceaa  in  de^  water.  Two 


I other  8itu«tl'’»ns  of  xaterast  are  that  of  a .ioastal  station  In  sI^How  water 

[ and  the  Intei-jaediate  situation  of  an  ocaao  station  in  relatively  shallow 

I wa.tv»r,  long  tern  records  exist  which  will  he  analysed  for  these  situations 

I I « in  the  mms  soanner  as  discussed  in  this  raeBaui  ecduifio  There  are  ^5  years  of 

I'  data  taken  at  a coasUJL  station,  namely  Scripps  pier  in  La  Jolla,  California, 

i Tiiere  are  soase  20  years  of  data  taken  at  each  of  two  stations  on  Queen 

Charlotte  Islruid,  lASTURA  and  ST.-  JAMBS,  which  because  of  their  expoeed 
iiature  car.  be  eonsidared  as  ocean  stations  in  shallow  water. 

' As  is  always  the  case,  it  would  be  desirable  to  obtain  wuch  longer  tine 

aeries  to  use  in  invaB’^lgatioos  for  trend,  end  for  Eore  extensive  auto- 
[ correl&tiaa  type  analyses.  The  long  records  aesticned  above  should  help 

1:  solve  this  prebiem. 


